Unstable particles as open quantum systems 
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We present the probability-preserving description of the decaying particle within the framework of quantum 
mechanics of open systems, taking into account the superselection rule prohibiting the superposition of the 
particle and vacuum. In our approach the evolution of the system is given by a family of completely positive 
trace-preserving maps forming a one-parameter dynamical semigroup. We give the Kraus representation for 
the general evolution of such systems, which allows one to write the evolution for systems with two or more 
particles. Moreover, we show that the decay of the particle can be regarded as a Markov process by finding 
explicitly the master equation in the Lindblad form. We also show that there are remarkable restrictions on the 
possible strength of decoherence. 

PACS numbers: 03.65.Yz 



I. INTRODUCTION 



Recently, the tests of Bell inequalities 1 1 ] in the system of 
correlated neutral kaons |2, 3] or B mesons [4] has attracted 
some attention. The crucial point in studying correlations in 
this system are the oscillations of the strangeness and bottom, 
respectively. However, the instability of kaons makes the anal- 
ysis of correlation experiments difficult. The state of the com- 
plete system is a superposition (or a mixture) of the states of 
the decaying particles and the decay products. The whole sys- 
tem undergoes a unitary evolution described usually in terms 
of quantum field theory. On the other hand, in correlation ex- 
periments of Einstein-Podolsky-Rosen-Bohm type lEll^], it 
is more useful to neglect the evolution of decay products and 
consider solely the decaying particles. Unfortunately, such 
a description within the framework of quantum mechanics re- 
ferring to the case with finite degrees of freedom leads to some 
difficulties. This is usually done by means of the Weisskopf- 
Wigner approach oUlL where the probability of detecting the 
particle is not conserved during the time evolution and, there- 
fore, the Hamiltonian in such theories must be non-Hermitian. 
Moreover, this formalism leads to some ambiguities when ap- 
plied to the description of correlation experiments. The reason 
is the probability loss caused by the decrease of the trace of the 
reduced density operator. This prevents one to calculate un- 
ambiguously the probability of finding the system in a given 
state after the projective measurement. Therefore, in our opin- 
ion, we need an approach enabling a description of the system 
that can be in two-particle states as well as the one-particle 
and even zero-particle states, which can arise during the time 
evolution (decay) of the initial system. It seems to us that the 
mentioned issues can be resolved by an assumption that the 
decaying particle can be found in a particle state as well as in 
the state of the absence of the particle, i.e., in the vacuum state 
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(it is not a vacuum in the sense used in quantum field theory, 
but rather in a sense used in |9[). 

In this paper, we give the probability preserving descrip- 
tion of the decaying particle within the framework of quan- 
tum mechanics of open systems. This approach is introduced 
in Sec. [II] where the evolution of the system is given by a 
family of completely positive trace preserving maps forming 
a one-parameter dynamical semigroup lUOt M U fl2i1 . Thus, 
in our approach the Hamiltonian is Hermitian and therefore 
the reduced density operator has a unit trace. We also find 
the operator-sum representation (Kraus representation) for the 
evolution of such systems, which immediately allows one to 
write the evolution for systems of two or more noninteract- 
ing particles. This is useful if we study quantum correlations 
between unstable particles. We would like to point out that 
we use the dynamical semigroup approach for the entire evo- 
lution of the unstable particle, not only for the description of 
its decoherence, as was done in 0[l3[ll[ll[li. Finally, 
in Sec.|ni|we study the restrictions on the possible strength 
of decoherence that arise as a side effect of completely posi- 
tive evolution of the system as well as we estimate the upper 
bound for the decoherence strength for K° and B° mesons. 



II. THE TIME EVOLUTION OF UNSTABLE PARTICLES 

In this section we discuss the evolution of unstable parti- 
cles, neglecting their spatial degrees of freedom. In order to 
sketch our approach, we begin with the discussion of the case 
of the neutral pion. This example is rather elementary, but 
it helps us to illustrate the main idea of our approach. Next, 
we go to K° (B°) mesons. The evolution of these particles is 
more complicated because of the phenomenon of transmuta- 
tion between K° and K° (B° and B~°). We shall regard them as 
open systems, it means that their evolution is not unitary, but 
it must be treated as a one-parameter family of quantum op- 
erations forming a dynamical semigroup. Consequently, the 
density operator of the system must obey the master equation 
rather than the von Neumann equation. 
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A. Unstable (pseudo)scalar particle 

In this section we describe briefly the evolution of an unsta- 
ble (pseudo)scalar particle, 7T°, and we show that this evolu- 
tion can be regarded as a family of amplitude damping quan- 
tum operations 11711 . 

The key point of the presented approach is that the system 
under consideration can be regarded as a two-level system: 
one can find the system in the particle state or in the vacuum 
state. Of course, this system must be an open one and we treat 
the decay products as a part of the environment. 

The space of states of the system is a direct sum of the 
Hilbert space of the particle H^o, spanned by the vector 
and the Hilbert space of the vacuum Ho, spanned by the vec- 
tor |0); i.e., H = H„o © Hq. Since we are looking for a de- 
scription of a decaying particle in terms of quantum mechan- 
ics with finite degrees of freedom (not in the language of field 
theory; cf. 11811 ') we assume that the decay process is a Markov 
process and it is not described in the dynamical manner (i.e., 
it is not governed by a Hamiltonian). 

We represent the vectors and |0) by 



10) 



(1) 



The time evolution of the system can be represented by the 
continuous one-parameter family of linear superoperators S t 
such that 



p(0=5,p(0), 



(2) 



where p(f) is the density operator of the system at the time 
t . These superoperators must be trace-preserving completely 
positive maps and they must form a one-parameter semigroup 
(see, e.g., 1 1 9l 12011 and refereces therein), i.e., 



S, 



te[<S f /3(0)]=tr[/3(0)] = l, 
SuS, 2 , Vfi,f 2 >0, 



(3a) 
(3b) 



and the map t i— > St is continuous in strong topology. 
We study the time evolution of the state of the system, 

p=pn\7t°) +p l2 \7t°) <0| + pf 2 |0) (jt°| + p 22 |0) (0|, (4) 

assuming that it is consistent with phenomenological 
Weisskopf-Wigner evolution \%\$, 



l*°(0> 



= g -?(im+r/2)| 7r 



(5) 



where m is the n° mass and F is its decay width. 

From l|5}, it follows that pu(t) = e~' r pn(0) and there- 
fore p 22 (f) = 1 — e ~' r Pn(0). Taking into account the linear- 
ity of S t , we can write Pn{t) as the time-dependent linear 
combinations of all the elements of the initial density ma- 
trix, i.e. pn{t) — j=iAij(t)pij(0), with the initial condi- 
tionsAi 2 (0) = 1 and all remaining A 's vanish at t = 0. 



Therefore the action of the map S t can be written as fol- 
lows: 



-tr 



S,p(0)=p(t) = 



Pn(0) 



i,j=\ 



. I a*.(Op*-(o) P22 ( f ) 

vj'=i 



(')py(o)^ 
/ 



(6) 



where 



P22W =p22(0)- 



Pn(0). 



To find conditions under which the map is completely 
positive, we use the Choi's theorem I19ll2lll . which states that 
S t is completely positive iff the corresponding Choi's matrix, 



Choi S, 



( e-' r A u (t) A u (t)\ 

A* n (t) l-e-T A* 21 (t) 

A 21 (f) A 22 {t 

\A* n (t) A* 22 (t) 1 / 



(7) 



is positive. This implies that 



|A 12 (0| 2 <^' r , 

|A„(0| 2 <(l-e- fr ) ( e -' r -|A 12 (r)| 2 ) 

A 21 (0=0, 
A 22 (f)=0. 



(8a) 
(8b) 
(8c) 
(8d) 



One can check by straightforward calculation that the com- 
position law ( 13^ leads to the conditions 

An(h+t 2 ) =A 12 (t2)A 12 (h), (9a) 
A u (*i +h) =A ll (t 2 )e- ,ir +A 12 (t 2 )A n (t l ). (9b) 

It is easy to see that the only possible solution of d9al fulfill- 
ing JHat and the initial conditions is 



Ai 2 (/)=^ r[(r+A)/2+i ' 11 , 



(10) 



where X > and jj, > 0. Taking into account the fact that the 
one-parameter semigroup must be Abelian, we get from J9bl >. 
and (JTHJ " 



A n (t) = 



z\e-' T -e 



zte 



-tr 



/[(r+X)/2+i M ] J ( when X^Torn^O, 
when A = r and ji = 0, 



(11) 



where z € C are such that the inequality J8bb is satisfied (for 
A = we have to put z = 0). Therefore, the most general form 
of the time-dependent density matrix is given by 
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Pit) 



-t[(T+X)/2-im] 



- tr Pn(0) 

p; 2 (P)+A*n(t)Pu(0) 



-t[(T+X)/2+im] 



pn(0)+A n (t)pu(0) 
Pn 



l- e -' r Pn(0) 



(12) 



r 



where the consistency with (0 requires fi — m. The parameter 
A is interpreted as the decoherence parameter. 

Since the evolution of the system fl!2i is given by a com- 
pletely positive map, it can also be written in the operator-sum 
form El 



p(t) = £E i (t)p(0)EHt), 

i=0 



(13) 



where the Kraus operators Ei(t) satisfy the condition 
J^LqEJ (t)Ei(t) = I. One can easily check that the Kraus op- 
erators leading to the evolution (11 2t are given by 

E (t)^e-^ r+x ^ 2+im ^\7t Q )(7Z°\ + \0)(0\, (14a) 

m = ^i-*- tr -^fif io><«°u d4b) 

= |0)(7r' 



E 2 {t)=e- ,T l 2 ^l-e-< x \% ){% 



,—tX 



(14c) 



From the operator-sum representation, using standard pro- 
cedures I23l 12411 . we can easily find the local form of the time 
evolution — the master equation in the Lindblad form lfl2ll . 

^ = -i[H,p(t)} + {K,p(t)} + ^L,p(t)Ll (15) 

where H is the Hamiltonian of the system, the operators L, 
are the Lindblad operators, and K — ~jYli=\E]U. For the 
density operator (II 21 . the Hamiltonian is 

H = m\%°){%X (16a) 
and the Lindblad operators are of the form 



L! = \/r(l-a)|0)<7r | 



(16b) 



L 2 = VI\% q ){k {) \ + ^=\Q){k\ (16c) 



where 



a = 



z| 2 [4;7t 2 + (r-A) 2 ]/(4rA), A ^Torm^O, 

A = r and m = 0, 

(17a) 



kl 2 /r 2 , 



\z\im+(T-Xy2\, A^rorm^O, ^ 
1 z, A = r and m = 0. 

Now, we take into account the fact that superpositions of 
the particle and vacuum are not observed in the nature. Con- 
sequently, there is no physical observable with nonvanish- 
ing matrix elements between vacuum state and particle state, 



which leads to the superselection rule. Therefore, the ele- 
ment P12 of the density operator @ does not contribute to 
the expectation value of any observable, and we can assume 
that Pn(t) — for any time l25ll . which implies that z = 0. 
Therefore in this case the decoherence parameter A becomes 
irrelevant, so we are free to put A = 0. Thus the density matrix 
describing the evolution of n° is 



Pit) 



'-' r pn(0) 

l- e - ,r p u (0) 



(18) 



The corresponding Kraus operators have the following form: 

Eoit) = e - f [ r / 2+, "V)(7r°| + |0)(0|, (19a) 



£ 1 (f) = V / l-^' r |0)^°|, 
while the generators of the master equation are 

6 = m\K°)(K°\, 
L!=\/T|0)<7r |. 



(19b) 

(20a) 
(20b) 



Note that the evolution of n° is thus simply the ampli- 
tude damping quantum operation JTvtl with the probability of 
damping depending on time, namely p = 1 — e~' r . If the ini- 
tial state of the system is p(0) = | tt°) (tt° | , then, as expected, 
the probability of detecting 71° at the time t is 

p(n ) = nip(t)\n°)(n Q \} = e-' r , (21) 

i.e., it is given by the Geiger-Nutall law. 



B. The time evolution of K° (5°) 



Now, we consider the case of a K° (B°) meson. This parti- 
cle needs special treatment because during the time evolution 
it transmutes into its antiparticle. Because both K° and B° 
mesons evolve according to the same scheme, hereafter we 
shall deal with K°, but the results are also valid for B° after 
appropriate changes of notation. 

The Hilbert space of the kaon-vacuum system fi K o © Ho 
is spanned by orthonormal vectors \K°), \K°), and |0), which 
are the eigenstates of the strangeness operator S: 

S\K°} = \K°), S\K°)=-\K°), S\0)=0. (22) 



These states (except of |0)) are not eigenstates of the opera- 
tor CP, where P is the space reflection and C is the charge 



conjugation. The CP eigenstates and \KS), 



CP\K 



0\ 



IK 1 



l 1 ). 



CP\K° 2 ) 



K$), CP\Q) = \0), (23) 
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are related to S eigenstates by 



where 



\K' 



1 



v/2 



(\K°) + \g°)) 



\K§) = -L (|*0) _ 1*0)) 



(24a) 
(24b) 



On the other hand, the time evolution operator is not diag- 
onal in the basis (I23> due to CP violation, but it is diagonal in 
the basis (see, e.g., 1201 1 



1 



l 



(\K°) + e\K°)), (25a) 



e|^> + l*2°» , (25b) 



with e being the complex CP-violation parameter, |e| « 
2.284 x 10 -3 - The time evolution in this basis is assumed 
to follow the Weisskopff-Wigner phenomenological prescrip- 
tion, 



\K$(t)) =e- t( - ims+rs ^\K°), 
|*0(0) =e"' (imi+ri/2) |«L), 



(26a) 
(26b) 



where F$ and Fl are decay widths of *°. and *°, respectively; 
m$ and mi are some parameters — their physical meaning is 
provided by the formulas ( I39l >. We would like to point out 
that these masses cannot be the eigenvalues of a Hermitian 
Hamiltonian because of a CP violation. Indeed, the basis ( 125 \ 
is no longer orthonormal, since 



2S(e) 
1 + lel 



8 L « 3.27 x 10" 



(27) 



and these states cannot be the eigenstates of a Hermitian op- 
erator. 

The most convenient way of analyzing the evolution of the 
density operator is to decompose it as follows: 

+ p* SL (t) \K° L ) (K$ I + p LL (t)\K° L ) (* L ° I + ~p L0 (t) \K° L ) (0| 
+ P5o(0|0> (K$ I + p£ (r) |0) \K° L ) + p 00 (0|0) (0| . (28) 

The superselection rule for a *° meson allows one to put 
Pso(t) — Plo(0 = for physical states. 

Because the basis d25l > is nonorthogonal, the matrix p(f) 
built from the coefficients of the decomposition d28l is not 
formed from matrix elements of the density operator p(f) in 
this basis, therefore one should be careful while operating on 
p(f), especially tr[p(f)] = 1 implies that 



tr\p(t)] = l-28 L $l[p SL (t)}. 



(29) 



(see AppendixlEl. 

Let us denote by p (t) the matrix formed by matrix elements 
of p(f) in the orthonormal basis {|*°), |*°), |0)}. From (|25} 
it follows that these two matrix representations of p(t) are 
connected by 



p(t) = Vp(t)V\ 



(30) 



V 



(l + lel 2 )- 1 ^ 
e(l + |e| 2 )- 1 /2 





e(l- 
(1- 



; |2)-l/2 
|2W/2 



(31) 







Now, we find the time evolution of the density operator in 
terms of the matrix p. Using d30i we can write 

p(f ) = S t p (0) = S t Vp (0) V f = VS t p (0)V f , (32) 



where p(f) = 5 r p(0) [cf. {30)]. 

The maps S t must be completely positive and must form a 
one-parameter semigroup. By virtue of (1321 the maps S t are 
composed from S t and the map ( 130b . One can easily check 
using Choi's theorem that the map d30t is completely positive, 
so the maps S t are completely positive iff the maps S t are also 
completely positive. It is much easier to find the conditions 
under which the latter maps are completely positive. 

After checking the conditions for complete positivity (see 
Appendix lAl and taking into account the superselection rule, 
we get the following time evolution of the matrix p(f): 

/ e-'^p ss (0) e -'( r +^-' A '")p 5L (0) 
p(t) = l e -t<r+x+*m)p* L ( ) e -'^p LL (0) 











PooWy 
(33) 

where Am = mi — ms (see later), T — (Tg + Tl)/2, and be- 
cause of ( I29> . 



Poo(0 



= (1- 
-25J 



-^)p ss (Q) + (l~e-'^)p LL (Q) 

^1 -^ (m - iAm) )p5L(0)]+Poo(0). 01) 



Note that, contrary to the case of % , the decoherence param- 
eter A is no longer irrelevant for evolution of physical states. 

The condition that the evolution of the density operator 
should be completely positive requires that the inequality 
[see jA9ll. 

8l(l-2e- rlr+ Vcos(tAm)+e- 2tir+ V)<(l-e- rrs ){l-e-' rL ), 

(34) 

must be valid for any t > 0. This inequality restricts the range 
of the parameters X, Tg, Tl, Am, and 8l- For the physical 
values of F$, Ti, Am, and Sl for* mesons and corresponding 
parameters for B° mesons, inequality d34t implies an upper 
bound on X (see Sec.llIU. 

We can write the evolution of the density operator given 
by d28i and d33i in the form of the operator-sum repre- 
sentation J 1 31 with the following Kraus operators (see Ap- 
pendixlBt: 



Eo(t) = 



1 



? -'[( r ^+^/ 2 +™d|*0)(*0| 



+e -,[(T L+ X)/2 + , mL ] lK } (A 0| _ §L r e -,[(T s+ X)/2 + im s ] |^ ) {R 1 



-?[(rz,+A)/2+i mL ]|^0w^0| 



-a/2 



|0> <0|, 



(35a) 



5 



E l (t) = 



I _ e -«(r+A-i'Am) I 



-&|0><je£|) ; 



(|0)<* s °| 
(35b) 



E 2 (t) 



' , 1 _ „-*(r+A--!Am)\ 

vT^^-4 2 . ... I |0)(*£| 



I _ e -t{T+X-i&mY 



-4 ( v^-^r,— _ • ) | 0)(A O| 



(35c) 



£ 3 (f) 



£ 4 (f) 



,—tX 



*°)(*°|-4|* S )(* L °|) 



(35d) 



,-tX 



(\K° L )(K L \-8 L \K L )(K° S \) 



^ 5 (0 = Vl-«- rt |0)(0|. 



(35e) 
(35f) 

Note that ( I34> ensures the reality of the square root in J35bl . 

The density operator p(t) fulfills the master equation dl5> 
with the following Lindblad operators: 



Jr s -8l\r+X-iAm\ 2 /r L 
i. =^ —(10X^1- 4|0)<*°|) 



r L -Sftr+x-iAm)/y/r£. 



'-\0)(K»\ 

/rZ-(r+A-/Am)/VrZ, 



1~4 2 



:|0)(*! 



(36a) 



(36b) 



u = 



(|*?)(* s |-4|*£)(*2|), (36c) 



1 _ 32 (l«?X«El-4l«2>(is5l) 



L 5 = VA|0)(0|, 



(36d) 



(36e) 



and with the Hamiltonian of the form (this is exactly the Her- 
mitian part of the Weisskopf-Wigner Hamiltonian): 



H = {'ns\K° s )(K Q s \+m L \K° L )(K\ 



-S L [(m - i AT/4) |J® <* L °| + (» + 1 AT/4) |*2) <*°|] } , 

(37) 

where Ar = F$ — Tl and m = (ra^, + ms)/2 is the mean *° 
mass (measured experimentally; m = 497.648 Me V/c 2 for 
*°, m B o = 5279.4 MeV/c 2 for 5° lEl l: 



m^o = (*°|H|* U ) =m, 
m^o = (*°|//|*°) = m. 



(38a) 
(38b) 



Note that m K o = mgo, as it is required by CPr theorem. Am is 
measured by an observation of *° flavor oscillation. This fi- 
nally gives us the interpretation of ms and m^, which appeared 
in J33i . as the expectation values of Hamiltonian in |*£) and 
|*°) states: 

m L = (*£ \H\Kl) = m + Am/2, (39a) 
m s = (Kl \H\K$) = m - Am/2. (39b) 

We would like to stress that the \K®) and |*£) are not eigen- 
states of the Hamiltonian d37l >. 

From ( I33> . it follows that the probabilities of detecting *° 
and *° are given by 



^o(f)=tr[p( f )|* )(*° 
+23? 



1 + 4 
2 

f(r+A-iAm) = 



'PSL 



{ e -^p ss (0)+ e -' r ^p LL (0) 
(40a) 



(0)]}, 



^o(0=tr[p( f )|*°)(* ' 



1 -^{e-' rs p SS (0)+e-'^p LL (0) 



-29? 



<• ' ' Psi (0)] } 



-r(r+A— /Am) = 



If the initial state is |*°), these probabilities are 
1 



Pgo(t) = 
pgo(t) = 



e- ,Ts +e- ,TL +2e- ,{T+X) cos{tAm) 

1 1 - 4 r 



(40b) 



(41a) 



4 1 + 4 



TL_ 2e -' (r+;i )cos(fAw) . 

(41b) 



The strangeness operator for *° is S = |*°)(*°| — 
|*°) (*°|, and its average is 



(S) = tv[ P (t)S] = 4 [e- ,Ts Pss(0) +e- trL p LL (0)} 



-25R 



-f(r+A-i'Am) = 



Psl(0) 



(42) 



Regardless of the initial state ]im t ^ ca (S) — 0, as expected, 
because in the limit t — ► oo we have vacuum only. If the initial 
state is |*°), we get 

1 ' .-t{T+X\ 



(S) 



1 



cos (f Am) + -^{e~ tTs 



.-*r, 



(43) 

Finally, we would like to point out that we are dealing 
with the Hermitian Hamiltonian and unit-trace density oper- 
ator. This assures us that there is no ambiguity in calculat- 
ing either conditional or joint probabilities for the results of 
measurements performed on the system. Moreover, we can 
unambiguously determine the states after the projective mea- 
surement using the standard quantum mechanical procedures 
(i.e., the von Neumann postulate of the state reduction). 



III. CP VIOLATION AND DECOHERENCE 

Now, let us analyze the inequality ( 134-t in more detail. Tak- 
ing into account that A > 0, we can treat J34l > as a quadratic 
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FIG. 1 : The discriminant of inequality 1341 as a function of time for 
K° and B°; in both cases A > 0, implying the validity of <44L 



inequality in e '( r +^). This inequality has real solutions pro- 
vided that its discriminant A fulfills 

A/8l = [(1 -e-' rs )(l -e- ,TL ) - 8l sin 2 (f Am)] > (44) 

for any t > 0, 
0.5292 x 10 10 s-\ 



where, according to 1 27]: Am — 
z s = l/r s = 0.8953 x 10- 10 ' 



I I / = 5.18 x l(T 8 s. Am B o = 0.502 x lO^s" 1 (Am B o w 
m B o -m B o), T = 1/r = 1.536 x 10~ 12 s, and 9?(%>)/(l + 

|%>| 2 ) = 0.5 x 10~ 3 . Fortunately, the inequality i44l holds 
for any f > for both K° and B° mesons because the first 
term rapidly grows from to 1 while the second one oscil- 
lates between and 5 2 ~ 10~ 6 (see Fig.^. Indeed, the series 
expansion of J44I is 



A/5 2 : 



[r 5 r L -5^(Am) 2 ]f 2 + o(f 



Therefore the condition 



5r< 



Am 



(45) 



(46) 



is necessary for the existence of the solutions of inequal- 
ity ( 134k This condition is satisfied for both K° and B° mesons. 

Now, let us analyze the restrictions imposed by the in- 
equality J34i on the range of the decoherence parameter A. 
From J34> we have, for any t > 0, 



cos (f Am) 



Va 
si 



< , 



-t{T+X) 



< cos (t Am) 



Va 



(47) 



FIG. 2: The allowed region for the decoherence parameter (gray) for 
K° and B° mesons: dashed curves are the upper and lower bound 
from I48i : the lower bound is always negative. The values of t+ and 
Amax for K° and 6° are given in the text. 



The left inequality gives a restriction only when its left- 
hand side is positive, which holds for K° and B° only for 



t < t+, where t+ w 7.18517 x 10 
1.53677 x 10" 15 s forfi . Therefore 
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s for K u and t+ 



A < — ln[cos(f Aw) - VA] -T, (48a) 

A > — \a[cos{tAm)+VK] -T, (48b) 

where the upper inequality must hold for < t < t + and the 
lower one for any t > 0; moreover, as mentioned earlier, A > 
0. These bounds are presented in Fig.[2]for K° and B°. Thus, 
for K° and 5°, 



< X < A ma 



inf 

0<f<r + 

1 



1 



In [cos (f Am) — VA] — T 



(49) 



r s r L 



<5 2 (Am) 2 - 



where the last equality comes from the first-order expan- 
sion. For K° we get A^ = 1.3629 x 10 11 s" 1 and for B° 
we get Amax = 6.5039 x 10 14 s~'. The experimental value 
for the decoherence parameter in an entangled K°K° system 



^3.80 . 



-3.30 



io^- 1 iEi 



is A = (1.84+|^) x 10~ 12 MeV = 2.80^ 
and for a B°B° system is A = (-47 ±76) x 10" 12 MeV 
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(-0.71 ± 1.15) x lO 1 ^- 1 jH and they fit in the allowed 
range. 

IV. CONCLUSIONS 

In this paper, we have found the probability-preserving de- 
scription of the decaying particle within the framework of 
quantum mechanics of open systems, taking into account the 
superselection rule prohibiting the superposition of the parti- 
cle and vacuum. It has been shown that some limitations of 
the Weisskopf-Wigner approach can be removed if we assume 
that the particle can be found in one of its possible states as 
well as in the state of the absence of the particle, i.e., in the 
vacuum state. In our approach the evolution of the system 
is given by a family of completely positive trace-preserving 
maps forming a one-parameter dynamical semigroup; thus the 
Hamiltonian is Hermitian and therefore the reduced density 
operator has a unit trace. It should be noted that we have 
used the dynamical semigroup approach for the description 
of the entire unstable particle, not only for the description 
of its decoherence, as in I3l 1 1 3l Il4l Il5l llfll . The advantage 
of the introduced approach is that there is no ambiguity in 
calculating either conditional or joint probabilities for the re- 
sults of measurements performed on the system. Furthermore, 
we can unambiguously determine the states after the measure- 
ment using the standard quantum mechanical procedures (i.e., 
the von Neumann postulate on the state reduction). We have 
also shown that there are restrictions on the possible strength 
of decoherence that arise as a remarkable side effect of com- 
pletely positive evolution of the particle. Moreover, we have 
found the operator sum representation (Kraus representation) 



for the general evolution of such systems, which allows one 
to write the evolution for systems with two or more particles. 
This is extremely useful if we study quantum correlations be- 
tween unstable particles. Moreover, we have shown that the 
decay of the particle can be regarded as a Markov process by 
finding explicitly the Lindblad form of the master equation for 
such a system. 
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APPENDIX A: TIME DEPENDENCE OF DENSITY 
MATRICES FOR K° 

In this appendix we show that the evolution given by 02} is 
the most general completely positive trace-preserving linear 
map that possesses the semigroup property and leads to the 
Weisskopff-Wigner evolution (I26> . This implies that we have 
well-defined p ss (t) = p ss (0)e- tT s and p LL {t) = p LL (0) e - ,r ^. 
Because kaons carry some quantum numbers like strangeness, 
we impose the superselection rule from the very beginning. 

The most general form of the evolution can be written in 
the form 



; Ass(0) 



S,p(0) 



I Aij(t)Pij(0) I B ij (t)p lj (t)\ 

iJ=S,L,0 iJ=S,Lfi 

. .£ A yW0y(o) e-^p^O) _ L L Qj(t)Pij(0) 

,' J £ ' »fj(0#S(0) I C*.(f)p*(0) hJ £ ' Dij(t)frj(p) 

\i,j=S,L,0 i,j=S,Lfl iJ=S,L,0 



(Al) 



The superselection rule causes that the only nonvanishing 
fi's and C's are Bjo(t), fio/(f)> Qo(0> an d Q)/(0> moreover, 
since p 00 (0 must be real, D u (f) = D* L {t), D os (t) = D* so (t), 
A)l(?) = Dfr.(t), and the other D's are real functions. The 
initial conditions are A^O) = 5yo(0) = Q,o(0) = 1 and the 
other functions vanish at t = 0; moreover, £>oo(0 = 1 because 



the vacuum must be a fixed point of this dynamics, i.e., 

p(0) = |0)(0|=*Vf >0: p(f) = |0)(0|. (A2) 
The corresponding Choi's matrix is 
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Ass(t) 








a 5l (0 








A so (0 


Bso(0\ 


















A 5s(0 





Cso(0 








D SS (t) 












Q(0 


D S o(t) 





















A L o(0 


Bw{t) 


a; l (0 










e -' r ^ 





aSl(0 





Qo(0 








%w 



















Aos(t) 


Bos(t) 





A 05 (f) 







A 00 (0 


fioo(0 







Cos(0 









A5o(0 





Coo(0 


V&(0 


q (0 




*2o(0 






*oo(0 


Q>(0 


1 / 



The positivity of this matrix requires that from the functions 
A's, fi's, and C's the only nonvanishing functions are AsL{t), 
Bgo(t), and Cuo(t). Moreover, the condition on the trace of 
p(f) implies that tr[p(f)] = tr[Vp(f)V" t '] = 1, so we have 

[e-' Ts +D ss (t)}p ss (0) + [e- tT t + D LL (t)}p LL (0) 
+ 2Vt{{8 I A SL (t)+D SL (t)}p SL (0)} 
+ 23i[D S o(t)p so (0)] + 25R[D i0 (f )Plo(0)] 

= Ass(0) +Pll(0) +25 L sR[p Si (0)]. (A4) 

Because this equation must be valid for any p (0), this implies 
that Djo(0 = for j 7^ 0, and 

£>oo(0 = l, (A5a) 

D ss (t) = l-e- tT s, (A5b) 

D LL (t) = l-e-' r \ (A5c) 

Offl(0 = &[1 -A S L(r)]. (A5d) 

Therefore, the matrix ( IA3I is positive, provided that 

|A 5L (?)| 2 <^' (r * +ri) , (A6a) 
|Bso(0| 2 <^' r S (A6b) 
|C L o(r)| 2 <^' r ^, (A6c) 

\B SQ (t)\ 2 e-'^ + \C L0 (t)\ 2 e-'^ + \A SL {t)\ 2 

- m[A SL (t)B* S0 (t)C w (t)] < e-'^ +r t\ (A6d) 



and 



\Dsdt)\ 2 <D ss {t)D LL {t). 



(A6e) 



Taking into account the composition law for the S t super- 
operator and ( IA6al i-( lA6cl . we get 



B so (t) = e-'^ 2+b +'^, 
C L0 (t)=e-'^ 2+c+i ^, 



(A7a) 
(A7b) 
(A7c) 



where a,fe,c,^i,^5,^£ > 0. Now, one can observe that the 
concrete time dependence of the functions Bso(t) and Cuo(t) is 
irrelevant, since the superselection rule requires that pso(0) = 



Plo(0) = and, therefore, pso(t) and Plo(0 still remain van- 
ishing (as required), so, without loss of any generality, we 
can choose a = b = c = A and Us = 111$, /Xl = mi, and then 
Hsl = Am. Thus, finally, we can put 



A SL {t)=e-'^ +T ^l 2+X - ihm \ 
B S0 (t) = e- t ^ 2+x+im s\ 
C i0 (f)=e-' [ri/2+A+imil , 



and, consequently, 



D SL (t) = 8 L (l- e -m+r L )/2 + i-iAm] 



(A8a) 
(A8b) 
(A8c) 



(A8d) 



Now, taking the solutions ( IA5a> -( lA5cl and ( IA8I . one can 
easily see that among conditions iA6\ for positivity of the 
matrix Choi 5, the only one that is not identically fulfilled 
is (IA6et: it can be written in the form 



§ 2 y_ e -t[(T s +T L )l2+X-iAm] 

for any t > 0. 



< (\-e- tTs )(\-e- tlL ) 



-tTr 



(A9) 



APPENDIX B: OPERATOR SUM REPRESENTATION FOR 



EVOLUTION OF K 







Let us denote 



1 8 L 0> 
V r V = I 8l 1 
1, 



We have 



tx\p{t)}=\x\p{t)g], 



(Bl) 



(B2) 



and therefore we can easily find Eq. \29\ . 

Now, let us find the operator sum representation for the evo- 
lution given by the map J33i . Let us define the set of matrices 
Ei(t), such that 



p(0 = ££K0p(0)£/(f 



(B3) 
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After a little algebra, we find that the matrices Ej(t) are 

/ e -t[{T s +k)j2+im s \ o o 

£«,(./ i = { e -t[{T L +X)l2+im L \ 

e-' x l 2 , 



(B4a) 



E 1 (t) = \ h-e-'^-Sl 



I _ e —t{T+X-iAm) 



2/0 0^ 



1 o y 
(B4b) 



£ 2 (f) = y/T^F^L 






1 _ e - f( r + A-iA, l Q 
6i l- e -'lL 1 




'1 0> 

^(fl^e-^Vl-r' 1 [ 

v Oy 

'o o o^ 



(B4c) 



(B4d) 



^(^e^'Vl-^' 1 10 

\0 Oy 

'0 0^ 



£ 5 (f) = \/l -e- tX [ 
.0 1 



(B4f) 



These matrices can be helpful in finding the Kraus operators 
Ei(t). Indeed, 

Y,V£i(t)p(0)£j{t)V^ =p(t) =Y,Ei(t) P (0)E}(t), (B5) 



where matrices on the right-hand side are written in the or- 
thonormal basis 123\ . This gives 



Ei{t) = VEi(t)V- 



(B6) 



More interesting than finding explicitly the matrices Ej(t )'s 
is finding the decomposition of Ei{t) into the sum 

E i (t)=£ i (t) ss \K s )(K s \ +E i (t) SL \K° s )(K° L \ +E,(t) SQ \K° s )(Q\ 

+£KOl S |^ )(^ I+£/(OllI^ >(^I+£/(OloI^°)(o| 

+ E i (t) os \0)(Kl\+E i (t) QL \0)(K L \+E i (t)oo\0)(0\. 

(B7) 



(B8) 



(B4e) Using Eqs. d25i and iB6\ . we get finally that 



E i (t)=V- 1 E i (t)V- lf =E i (t)g-K 
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